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■ High- resolution radar images from the Magellan spacecraft have revealed the first details of 

the morphology of the Lavinia Planitia region of Venus. A number of geologic units can be 
distinguished, including volcanic plains units with a range of ages. Transecting these plain, over 
much of the Lavinia region are two types of generally orthogonal features that we interpret to be 
compressions! wrinkle ndges and extensions! grooves. The dominant tectonic features of Lavinia 
are broad elevated belts of intense deformation that transect the plains with complex geometry, 
they are many tens to a few hundred kilometers wide, as much as 1000 1cm long, and elevated 
hundreds of meters above the surrounding plains. Two classes of deformation belts are seen in 
the Lavinia region. Ridge belts” are composed of parallel ridges, each a few hundred meters 
in elevation, that we interpret to be folds. Typical fold spacing* are 5-10 km. "Fracture belts” 
are dominated instead by intense faulting, with faults in some instances paired to form narrow 
grabens. There is also some evidence for modest amounts of horizontal shear distributed aooss 
both ridge and fracture belts. Crosscutting relationships among the belts show there to be a 
range m belt ages. In western Lavinia in particular, many ridge and fracture belts appear to bear 
a relationship to the much smaller wrinkle ridges and grooves on the surrounding plains: Ridge 
morphology tends to dominate belts that lie more nearly parallel to local plains wrinkle ridges, 
and fracture morphology tends to dominate belts that lie more nearly parallel to local r Ui ~ 
grooves. We use simple models to explore the formation of ridge and fracture belts. We show 
that convective motions in the mantle can couple to the crust to cause horizontal stresses of a 
magnitude sufficient to induce the formation of deformation belts like those observed in Lavinia. 
We also use the small-scale wavelengths of deformation observed within individual ridge belts to 
° Wer Umit <m th * vem “ iM1 dermal gradient in the Lavinia region at the 


Introduction 

Of the many linear tectonic features known to be 
present on Venus prior to the Magellan mission, most 
attention was devoted to the spectacular mountain 
belts found in the Maxwell Montes and elsewhere [e.g., 
Campbell et ai, 1983; Barsukov et al, 1986; Basilevsky, 
1986]. It has also been known, however, that there 
is considerable evidence for tectonism on the venusian 
plains. The “ridge belts” described on the basis of 
Arecibo and Venera 15/16 data are particularly good 
examples. In Venera 15/16 images they are most abun- 
dant in Atalanta Planitia, a lowland region in the north- 
ern hemisphere. Like the highland mountain belts, 
ridge belts were inferred to consist of subparallel ridges 
and grooves. However, they appeared to be elevated 
comparatively little above the surrounding plains. The 
ridges and grooves that form them are approximately 
parallel to the margins of the belts, and can be continu- 
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? us strike for hundreds of kilometers. Sp&cings of 
individual ridges, as determined from the low-resolution 
Arecibo and Venera 15/16 data, were seen to be typ- 
ically a few tens of kilometers. The belts in Atalanta 
Planitia are themselves crudely parallel, with typical 
belt-to-belt spacings of 300-400 km. In some locations, 
both in Atalanta and elsewhere, individual ridges and 
furrows not organized into belts also were observed. A 
compressions! origin for the ridge belts was the most 
widely accepted hypothesis prior to Magellan [e.g., Bar- 
sukov et ai, 1986; Crvmplcr et ai, 1986], although 
an extension al origin also was suggested \Sukhanov and 
Pronin, 1989]. 

The first deformation belts on the venusian plains 
observed by the Magellan spacecraft have been those 
on Lavinia Planitia, in Venus’s southern hemisphere. 
Solomon et al. [1991] presented a preliminary descrip- 
tion of some of the belts of Lavinia. They noted that 
there were two classes of belts and concluded that both 
had a compressions! origin. In this paper we discuss the 
belts of Lavinia in detail, using improved data and ex- 
panded imaging coverage. Because the Magellan images 
have shown the belts to have significantly greater mor- 
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phologic diversity than was apparent in previous data 
sets, we have followed Solomon et al. (1991] by adopting 
the general term “deformation belts” for these features, 
retaining the term “ridge belt” for one of the morpho- 
logic subtypes. 

In preparing this paper, we have primarily made use 
of three Magellan data products: full-resolution mo- 
saicked image data records (F-MIDRs), compressed mo- 
saicked image data records (Cl-MIDRs), and global to- 
pographic data records (GTDRs). F-MIDRs are syn- 
thetic aperture radar (SAR) image mosaics that retain 
the nill spatial resolution of the radar system, sampled 
at 75 m/pixel. Cl-MIDRs contain the SAR data re- 
sampled at 225 m/pixel. Altimetry is taken from GT- 
DRs, which contain measurements of local planetary ra- 
dius made by the Magellan altimeter, sampled to a uni- 
form planetary grid and spatially filtered for smoothing. 
The along-track spatial resolution of the altimeter varies 
from about ~8 km near the equator to ~20 km at high 
latitudes, and the across-track resolution varies simi- 
larly from —13 to ~30 km [Pettengill et al., 1991]. Due 
to a software error, the altimetric results for Lavinia and 


elsewhere presented by Solomon et al. [1991] are incor- 
rect at small spatial scales, and we present improved 
altimetric data here. 

Geologic Units 

Most of Lavinia Planitia is seen in the Magellan Cl- 
MIDR 45S350, shown in Figure 1. Several geologic 
units, defined on the basis of radar brightness, small- 
scale texture, characteristics and abundance of super- 
posed structural features, and apparent relative ages, 
can be distinguished in this area. Because the SAR inci-’ 
dence angle in this region averages about 25®, the radar 
return relates primarily to wavelength-scale roughness 
and, less commonly, to meter-scale and larger slopes; 
the effects of compositional variations should be minor! 
Consequently, many, and perhaps most, of the geologi- 
cal units are distinguished by differences in surface char- 
acteristics rather than differences in inherent material 
properties. For some units, such as lobate flows, radar 
units very likely correspond to those defined using more 
conventional geological criteria. Other units, however, 
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may reprint similar materials that have experienced 
different structural evolutions, with consequently differ- 

en r^ teX J UrCS at cent * ine ter end meter scales. 

We divide the terrains of Lavinia Planitia into four 
broad groups: complex ridged terrains (CRT), textured 
terrains, regional plains, and digitate plains. We have 
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Fig. J. C©raplex ridged terrain, exposed through younger plains 
nutenak. Width of the image is 190 km. 
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defined one CRT unit, one textured unit, two regional 
plains units, and one digitate plains unit. The follow- 
ing paragraphs describe the units in approximate order 
from oldest to youngest. 

Complex Ridged Terrain 

CRT occurs in Lavinia Planitia as scattered inliers 
surrounded and embayed by younger plains deposits 
(Figure 2). These inliers are characterized by a strong 
linear trend of kilometer-scale or narrower ridges and 
troughs, with a weaker secondary trend of lineaments 
at a high angle to the strong primary trend present lo- 
cally. CRT is much brighter than surrounding plains 
on the SAR images unless the plains are intensely frac- 
tured. 

Textured Terrain 

Many areas in Lavinia are as radar-bright as CRT 
but do not exhibit the same kilometer-scale ridge-and- 
trough pattern. In places, it is possible to resolve a 
very regular fabric of bright lines at a scale of about 
400 m, and it is this fabric that gives this material its 
textured appearance (Figure 3). At the scale of the 
Magellan SAR images, this fabric can be considered 
penetrative. The structural cause of the fabric cannot 
be determined, but a uniformly spaced set of joints is 
one possibility. Textured terrain occurs in two distinct 
physiographic forms: (1) as plains that are significantly 
brighter than adjacent regional plains materials, and (2) 
as long, narrow ridges generally grouped together into 
ridge belts (see below). Crosscutting and embayment 
relationships in some instances indicate that textured 
terrain is younger than CRT and older than regional 
plains. Nevertheless, both textured terrain and regional 
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plains arc of variable brightness, and therefore there are 
places where it can be difficult to define the contact be- 
tween them. The penetrative fabric and, indeed, the 
overall textured aspect are not evident in places where 
textured terrain is intermediate in character between 
well-defined textured terrain and bright mottled plains. 
Regional Plains 

The areally dominant terrain type in Lavinia Planitia 
consists of moderately radar-dark to moderately radar- 
bright plains. Regional plains materials are most likely 
of volcanic origin, but morphologic forms resembling 
flows are rare. 

Mottled plains. These are moderately bright to locally 


bright plains with a “blotchy” or mottled texture at 
kilometer and smaller scales (Figure 4). They are gen- 
erally characterized by abundant superposed wrinkle 
ridges (see below). In places, mottled plains include 
large numbers of small shields or domes from a few to 
about 20 km in diameter. Locally, it is possible to de- 
fine two members of mottled plains that differ in radar 
brightness. At many such locations it is the brighter 
and apparently older of the two members that includes 
the small shields and domes, which are surrounded by 
the darker and apparently younger member (Figure 4). 

Dark plains. These are plains that are darker, have 
more uniform radar backscatter, and have fewer super- 
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posed wrinkle ridges than mottled plains. Dark plains 
commonly occur as small patches adjacent to or within 
fracture belts (see below), where they are superposed 
on the fractures. In other places, however, dark plains 
appear to grade into the dark member of the mottled 
plains (Figure 4). 

Digitate plains 

This unit consists of complexes of digitate to locally 
lobate flows, many attaining lengths of hundreds of kilo- 
meters (A in Figure 1). Most of these flows are brighter 
on the SAR image, and hence probably rougher at cen- 
“ " time ter scales, than all other plains units. Feeder chan- 
nels (or collapsed lava tubes) are apparent in several 
places, and many flow complexes can be traced ups- 
lope to sources of apparent volcanic origin. Superposed 
wrinkle ridges are sparse, and flow patterns clearly in- 
dicate that adjacent deformation belts already existed 
as topographic ridges when these flows were emplaced. 

Tectonic Features 

_ Wrinkle Ridges 

Most of the plains areas of Lavinia contain long, nar- 
rpwy sinuous linear features that appear brighter than 
the background plains (Figure 4). These features are 
generally less than 1 km wide and a few tens of kilome- 
ters long, although some reach widths of several kilome- 
ters and lengths in excess of 100 km. Typical spacings 
range from several to about 20 km and locally up to 


50 km. Where the topography of these features can be 
inferred from brightness variations in the SAR images, 
they are clearly seen to have positive relief, although 
most appear as simple sinuous bright lineaments. We 
interpret these features as compressions! ridges, anal- 
ogous to the wrinkle ridges commonly seen on the lu- 
nar maria and Martian ridged plains [e.g., Plescia and 
Golombeky 1986; Watters t 1988]. 

The wrinkle ridges of Lavinia Planitia show strong 
preferred orientations within subareas of the region 
(Figure 5«). A very strong NNE-SSW trend is apparent 
in the western portion of the map, while the trend is N- 
S in the central portion and NNW-SSE in the northeast 
portion. These trends are also evident in rose diagrams 
of ridge segment orientations (Figure 5t). 

Wrinkle ridges appear to have formed throughout 
the evolution of plains, because all plains units have 
them, including digitate plains. However, there is a 
clear decrease in the number of ridges from mottled 
plains to dark plains to digitate plains; that is, there 
is a decrease in ridge concentration with decreasing in- 
ferred age. Furthermore, the trends of the ridges ap- 
pear related to locality rather than stratigraphic age 
because the dominant trend of the sparse ridges on dig- 
itate plains is similar to the dominant trend of the more 
abundant ridges on adjacent regional plains. These re- 
lationships suggest that ridge formation is an ongoing 
process related to local stress fields that did not change 
significantly during the deposition of the various plains 
materials. Although local stress fields varied geograph- 




Fig. 5*. Map of wrinkle ridges (line* on plains) and deformation belts in Figure 1. For fracture belts, fracture 
trends are shown by dashed line patterns; for ridge belts, ridge trends are shown by lines with diamond symbols. 
Shaded regions are complex ridged terrain. 
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Grooves 

Over much of the western part of Lavinia Planitia are 
WNW trending radar-bright lineaments. They are geo- 
metrically distinct from the wrinkle ridges in that they 


are both longer and substantially straighter. In many 
instances they are too narrow to be resolved as anything 
but radar-bright lines, but where they can be resolved 
they are seen to be narrow grooves (Figure 6). Lengths 
from 25 to 75 km are common in some areas, while 
lengths of 150 to 200 km are common in others. A few 
grooves have lengths exceeding 500 km. Spacings are 
also variable. The longest grooves (lengths exceeding 
200 km) show typical spacings of 30 to 100 km. Heav- 
ily deformed regions show typical spacings of 1 to 3 km. 
Groove widths are generally about 300 to 700 m. There 
are no instances of these grooves cutting digitate flows, 
but it is not clear if this is because the flows are younger 
or simply because these grooves were never formed in 
the eastern portions of Lavinia where such flows occur. 
The grooves clearly are younger than mottled and dark 
plains. Abutting relationships with wrinkle ridges sug- 
gest that some ridges are older than adjacent grooves, 
some younger. We interpret these grooves to be narrow 
grabens. 

A map of grooves in Lavinia Planitia shows the par- 
allel nature of the grooves in the western part of the 
region (Figure 7). It is noteworthy that in all areas 
where grooves are present, their orientations are per- 
pendicular to those of the wrinkle ridges. Thus, a sig- 
nificant section of Lavinia Planitia is characterized by 
a distinctive pattern of orthogonal compressional and 
extensional features on the plains between the defor- 
mation belts. 



Fig. 6. Groove* on the plain*. Width of the image is 205 km. See Figure 1 for location. 
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Fig. Tb. Roae diagram, of groove orientation, in aix mbarea. of 
rigure 1. Normalization is as in Figure 5b. 

Ridge Belts 

The most obvious tectonic features of Lavinia Plani- 
tia are the prominent deformation belts that transect 
much of the surface. These belts take two very distinct 
forms which we call ridge belts and fracture belts. An 
example of a ridge belt is shown in Figure 8. The ridges 
that compose them can vary considerably in width both 
from one to the next in a given belt and along the 
length of a single ridge. The maximum width observed 


is about 10 km, and more typical widths are a few kilo- 
m *!f re - , Some Hdges appear symmetric in cross section, 
while others appear steeper on one flank than the other. 
A smooth, archlike proflle is most common, but some 
ndges show a narrow, rugged secondary ridge superim- 
posed on the crest or on one flank of the arch as is 
common for junar wrinkle ridges. The ridges are typi- 
cally rather sinuous and commonly bifurcate and merge 
along strike, producing a complex anastomosing pat- 
tern. borne ridges appear to consist of textured terrain 
material, but m most places the unit involved is not 
apparent. 

It would be desirable to construct detailed topo- 
graphic profiles across ridge belts, but the individual 
ndges are too narrow to be resolved by the Magellan 
altimeter, and the single SAR view obtained during 
the first mapping cycle does not permit stereo analy- 
sis. There is, however, one approach that can be used 
to construct such profiles. One of the belts (Figure 8) is 
®f vera ^ f^e tectonic grooves discussed above. 
While these grooves appear very straight where they 
he on the adjacent plains, they take on a distinctive 
cuspate curvature as they cross each ridge in the belt. 

It is likely that this curvature is an effect of the radar 
viewing geometry; the low radar incidence angle at this 
latitude (about 23*) results in the groove, as it passes 
across the elevated ridge, appearing displaced in the di- 
rectjon of the spacecraft. 

If we make the assumption that the grooves are ac- 
tually straight in map view as they cross the ridges, 
we can use the apparent local deviation of the gr<x>ve 
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trend, plus a knowledge of the radar viewing geome- 
try, to construct a topographic profile. Such a profile 
is shown in Figure 9. Most of the ridges in this profile 
are fairly symmetric, with archlike profiles, widths of 
several kilometers, and heights of 200-300 m. (This ap- 
proach to producing a topographic profiles is, of course, 
fundamentally limited by the correctness of the assump- 
tions on which it is based. The detailed heights and 
shapes of the ridges in the profile may be fairly accu- 
rate, since grooves elsewhere tend to be very straight 
on scales comparable to the ridge width. The overall 
height of the belt in the profile, however, should not 
be regarded as accurate, since subtle curvature of the 
grooves over length scales of many tens of kilometers is 
common.) 

The individual ridges within ridge belts differ little 
in appearance from the largest wrinkle ridges, except 
that they tend to be less sinuous. Because both types 
of ridges are inferred to result from compressions] de- 
formation, this similarity is not too surprising. Where 
ridge belts deforming textured terrain are embayed by 
plains materials with wrinkle ridges, large, young wrin- 
kle ridges can be distinguished from older ridges of the 
belt. The possibility must be entertained, however, that 
ridge belts can evolve from wrinkle ridges and thus be 
younger than the adjacent regional plains at some loca- 
tions. 


Planitia is the fracture belt (Figure 10). In contrast 
to ridge belts, fracture belts are dominated by a com- 
plex pattern of linear to arcuate faults and fractures. 
Many faults appear singly, but others are paired to 
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Fracture Belts 

The other class of deformation belts found in Lavinia 


Fig. 9. Inferred topographic profile across the ridge belt in Figure 
9, constructed from the apparent displacement of a groove as it 
crosses the belt. 
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form grooves with widths ranging from a few kilometers 
down to the resolution limit of the images. The faults 
commonly display complex anastomosing and crossing 
patterns that indicate repeated or progressive deformar 
tion. Where fracture belts change trend, or where belts 
bifurcate or merge, the pattern of faults is especially 
intricate. 

In general, fracture belts deform mottled plains ma- 
terials but not dark plains materials. The obvious di- 
version of flows and of channels in digitate flow fields 
around some fracture belts in eastern Lavinia (B in Fig- 
ure 1) and the truncation of the fractures by flows in- 
dicate quite clearly that these fracture belts are older 
than the digitate flow fields. 

A noteworthy aspect of a few of the fracture belts is 
that in some parts of them two distinct scales of defor- 
mation are observed. Faults are spaced very closely 
(typically a few hundred meters) almost everywhere 
they are present in fracture belts. However, in a few 
locales these closely spaced faults are concentrated in 
bands of intense deformation that are separated by 
nearly undeformed materials. The bands typically are 
spaced 20-30 km apart. It was reported by Solomon et 
al [1991] that these bands of intense deformation co- 


incided with local highs in the altimetry, and this was 
interpreted to mean that the faulting might be a conse- 
quence of brittle deformation across the crests of anti- 
clines. However, as mentioned above, this inference was 
based on incorrectly processed altimetric data. In Fig- 
ure 11 we present the same altimetric scan across one 
of the fracture belts, Hippolyta Linea, that was shown 
by Solomon et al [1991], here using the reprocessed 
altimetry. It now appears that the intensely fractured 
bands within the belts are not higher than the adjacent 
materials; if anything, they may be lower. Definitive 
treatment of this problem will have to await acquisition 
and processing of stereographic data, but at present the 
idea that the bands coincide with anticlinal crests ap- 
pears uncertain. An alternative possibility is that the 
bands are regular concentrations of extensional defor- 
mation caused by necking instabilities. 

Much of the faulting that is so prominent in the frac- 
ture belts is similar in appearance, though much greater 
in density, to the faulting on the plains that forms the 
grooves there. In fact, there are a number of instances 
(e.g., Figures 1 and 10 and Figure 17 below) where fault- 
ing in fracture belts continues onto the plains, at much 
lower density, to form widely spaced grooves. As is the 



Fig. 10. A fracture belt. Width of the image is 260 km. See Figure 1 for location. The altimeter ground track for 
orbit 415 is shown. Arrows correspond in position to arrows in Figure 11. 
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Fig. 11. Altimetry data from orbit 415. Arrows corre s pond in 
position to arrows in Figure 10. 

case for the ridge belts and the wrinkle ridges, then, 
there appears to be a close relationship between fault- 
ing in fracture belts and the formation of grooves. 



Altimetry 

Magellan altimetry data have also provided the first 
detailed look at the topography of Lavinia Planitia. 
Figure 12 presents a topographic contour map for the 
Lavinia region. Overall, Lavinia is a: broad, low plain, 
lying at a typical planetary radius of about 6050.5 km 
The regional topography rises gradually by more than 
1 km to the northwest, northeast, and southeast. The 
general rise in elevation to the east contbues well be- 
yond the edge of our map and culmbates b a Ibe of 
several large coronae south of Alpha Regio [Stofan et a/., 
this issue]. A prombent characteristic of the altimetry 
of Lavinia is that both types of deformation belts are 
elevated above the surrounding plains. 

The fracture belts appear b general to reach higher 
elevations than the ridge belts. A number of the frac- 
ture belts reach elevations that are as much as 1 km 
above the level of the adjacent plains, and elevations 
500 m above the plains within fracture belts are the 
norm. Particularly high elevations commonly are found 
in complex and intensely deformed regions present at 
bends and junctions in fracture belts. In contrast, max- 
imum relief of 200-300 m is more typical for the ridge 
belts, and some of the ridge belts are not evident at 
all b the altimetric data. Some care must be exer- 
cised b interpretation of the altimetry, however, as de- 
termination of surface elevation for a given altimeter 
pulse can depend sensitively on both the large-scale and 
small-scale surface roughness within the pulse footprint. 
Roughness on a large scale can cause multiple echoes 
that severely complicate the range-finding algorithm. 
Roughness on a small scale can result in a surface scat- 
tering law that deviates from the one assumed in fitting 
the shape of the leading edge of an echo, again affectbg 
the range-findbg process. In any case where two ter- 
rains of significantly different roughness are observed, 
systematic variations in elevation between them must 
be treated with caution. 


Lavinia: Contour Interval * 100m 
Latitude: 44.2S - 50.9S 
Longitude: 347.6E - 356.9E (at north edge) 

Fig. 12. Topographic contour map for the area covered by Figure 
1, in ainueoidal equal-area projection. Elevation, are referenced 
to a planetary radius of 6050.0 km. 


Geologic Relationships 

Evidence for Horizontal Shear in Deformation Belts 

Solomon et al. [1991] presented two pieces of evi- 
dence for modest amounts of horizontal shear parallel to 
the axes of deformation belts. These are shown b Fig- 
ures 8 and 10. In Figure 8, an old set of plains grooves 
exhibits an S-shaped bend consistent with distributed 
left-lateral shear as it crosses the ridge belt. In Fig- 
ure 10, the fracture belt exhibits a number of instances 
where paired faults form rhombohedral downdrops, b 
a manner similar to what is seen b terrestrial settings 
where extension is coupled with shear. In this instance, 
the bferred sense of horizontal shear across the belt is 
right-lateral. 

There are several other instances in Lavinia Plani- 
tia where patterns developed within deformation belts 
can be bterpreted as very large-scale analogs of famil- 
iar kinematic bdicators seen in terrestrial shear zones. 
One of these (Figure 13) shows a series of S-shaped mi- 
nor folds along the lengths of bdividual ridges within 
a ridge belt. S folds of this geometry are indicative 
of left-lateral (sinistral) movement and in this bstance 
would imply sinistral shearing parallel to the length of 
the ridge belt. The belt showing the S folds termi- 
nates to the SSW against a contact with mottled plabs, 
and there is no indication of sinistral slip within these 
plains. Consequently, the inferred belt-parallel sbistral 
slip must be older than the mottled plabs. It either rep- 
resents a late stage in the development of the ridge belt 
or implies that the ridge belt is transpressional (com- 
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Fig. 13. S-ihaped folds in a ridge belt, with orientations consis- 
tent with left-lateral shearing of the belt. Width of the image is 
130 km. See Figure 1 for location. 


pression plus sinistra] shear) rather than purely com- 
pression al in origin. 

A second example of a possible kinematic indica- 
tor (Figure 14) involves the distortion of well-developed 
penetrative fabric in textured terrain as the fabric im- 
pinges on a throughgoing, strongly lineated boundary of 
a ridge belt. This boundary abrutly truncates all struc- 
tures within the belt and is interpreted to be a shear 
zone. If this inference is correct, the distortion of the 
fabric at this location implies a right-lateral (dextral) 
sense of shear. 

Age Relationships Among Belts 

The deformation belts in Lavinia did not all form 
contemporaneously. Instead, they show some signifi- 
cant variability in age, both with respect to the plains 
materials surrounding them and with respect to one an- 
other. As noted above, a number of belts, particularly 
in western Lavinia, appear to deform the same mate- 
rial that constitutes the adjacent plains, with little or 
no subsequent volcanism. In the ridge belts, the tran- 
sition from dark or mottled plains adjacent to the belt 
to the textured terrain of the belt is not necessarily a 
stratigraphic one in all cases; we cannot rule out the 
possibility that in some instances the texturing simply 
resulted from the folding of the plains material. How- 
ever, there are other very clew instances where belts are 
embayed and partially buried by lavas that postdate the 
belt deformation. The N-S fracture belt at C in Figure 
1 is one of these, and it may have once been contiguous 
with another aligned and partially buried N-S belt at 

C. 

An example of a clear crosscutting relationship in- 
volving two belts is shown at D in Figure 1. Here a N-S 
belt has been largely buried by lavas, and these lavas 
have in turn been deformed to form an E-W belt that 
crosscuts the N-S one. This case in particular demon- 
strates that belt formation in Lavinia took place over an 
extended period of time, during which which belt for- 
mation episodes were interleaved with plains volcanic 
events. 

Regional Tectonic Trends 

As noted above, there are two classes of deformation 



Fig. 14. Penetrative fabric in textured terrain that is distorted 
with a geometry consistent with right-lateral shear across a ridge 
belt. Width o I the image is 125 km. See Figure 1 for location. 


belts in the Lavinia region, and each appears related 
to tectonic features on the intervening plains. Specifi- 
cally, ridges of the ridge belts are morphologically sim- 
ilar to and in some cases contiguous with plains wrin- 
kle ridges, and faults of the fracture belts are similar 
to and in some cases contiguous with plains grooves. 
In some parts of Lavinia, however, these relationships 
may be more fundamental, as suggested by Figures 5 
and 7. The morphologic form taken by each deforma- 
tion belt, ridge belt or fracture belt, appears related 
in many instances to the orientation of the belt with 
respect to the orthogonal tectonic features on the ad- 
jacent plains. Belts whose orientations are closest to 
those of the plains grooves in their vicinity tend to be 
fracture belts, while those whose orientations are closest 
to those of local wrinkle ridges tend to be ridge belts. 
Even more notably, there are instances, as in Figure 15, 
where a single belt changes from one trend to the other, 
and in so doing also changes tectonic style from one belt 
type to the other. In addition, there are some belts (e.g., 
Figure 16) where both ridges and fractures are common 
within the belt and where the ridges lie parallel to local 
plains wrinkle ridges and the fractures lie orthogonal to 
them and parallel to local plains grooves. 

There are exceptions to this relationship. The most 
obvious one is the N-S fracture belt at C in Figure 1. 
In this same belt, some ridges also are present and lie 
parallel to fractures, rather than perpendicular to them. 
Other belts departing from the general relationship oc- 
cur farther to the east. Despite these exceptions, the 
relationship described above appears sufficiently preva- 
lent, particularly in western Lavinia, that it must be 
accounted for in some fashion. The belts in central and 
eastern Lavinia that deviate from the general pattern 
are in most cases very clearly buried and embayed by 
younger volcanic deposits and perhaps are formed in a 
stress field distinct from the one responsible for most of 
the wrinkle ridges and grooves. 
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Discussion 

Geologic Evolution of Lavinia Planitia 

The description above of geologic units and tectonic 
features in Lavinia Planitia leads us to some preliminary 
interpretations: 

Regional tectonic pattern on the plains. The orthog- 
onal pattern of compressional wrinkle ridges and exten- 
sional Proves that transects the plains units over much 
of Lavinia Planitia is remarkable in its regularity and in 
both its spatial and temporal extent. These character- 
ises indicate that at least the near-surface materials 
of Lavinia have been subjected to a stress field that has 
been generally uniform, both spatially and temporally, 
for much of the region’s history. This conclusion is not 
unique to Lavinia, since a variety of regular tectonic 
patterns are observed on the plains elsewhere on the 
planet [Solomon et al, 1991, this issue). 

It is tempting to think of the wrinkle ridges and 
grooves on the plains as the simple manifestations of 
thrust faulting and normal faulting, respectively, and 
this interpretation may well be correct. If so, it has 
some interesting implications. A common approach to 
predicting faulting style from stress field descriptions 
is that of Anderson [1951]. With the Anderson fault- 
ing criteria, however, no single set of values of the three 
principal stresses <rj , < 72 , and <73 simultaneously predicts 
orthogonal grabens and thrust faults. Applying An- 
derson criteria rigorously, one can obtain either thrust 


faults or grabens from a given stress field orientation 
but not both. ’ 


orthogonal thrust faults and grabens may both be 
formed, according to Anderson criteria, if the two stress 
components in the plane of the surface maintain their 
relative values but vary appropriately with respect to 
the vertical stress. One way this could take place would 
be for the in-plane stresses to vary relative to the ver- 
tical stress over time. Another would be if fault initia- 
tion for the two classes of features took place at differ- 
ent depths [Golombek, 1985], near the surface for thrust 
faults and at depth for grabens, so that in the former 
case the vertical stress is <73 and in the latter case it is 

<?i. (The depth is limited, however, since many grabens 
are very narrow.) 

In either of these scenarios there is a stress state, 
intermediate in either time or depth, for which Ander- 
son criteria would predict strike-slip faulting. Clear ev- 
idence for strike-slip faulting is rare on Venus (and in- 
deed on most of the terrestrial planets save Earth), and 
is lacking entirely on the plains of Lavinia. One expla- 
nation offered for the absence of strike-slip motion on 
other planets is that layered mechanical discontinuities 
concentrate stress, and hence fault initiation, at certain 
depths (and hence at certain confining stresses) where 
faults with styles other than strike-slip tend to be ini- 
tiated [Golombek, 1985]. Another possibility is simply 
physical restriction of lateral movement. A planet has 
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Fig. 16. A deformation belt that exhibits both ridges and fractures, with approximately orthogonal orientations. 
Width of the image is 240 km. See Figure 1 for location. 


a free surface and a ductile region at depth, so thrust 
and normal faulting within the strong upper crust can 
always be accommodated by vertical motions, and the 
stress at the base of a fault cutting the entire strong 
upper layer can relax by ductile flow. On a one-plate 
planet or within the interior of a single lithospheric 
plate, however, there may be restrictions on the extent 
of horizontal strain, and stress relaxation at the ends of 
strike-slip faults will not be comparable to that for dip- 
slip faults. It may be that such restrictions inhibit most 
large-scale (i.e., easily observed) strike-slip motions in 
Lavinia. Schultz and Zuber [1991] have recently sug- 
gested that in situations where the lithosphere is ver- 
tically loaded, tensile failure rather than shear failure 
may relieve much of the stress and may suppress the 
zone in which shear failure criteria predict strike-slip 
faulting. 


Deformation belts . In contrast to pre-Magellan ex- 
pectations, we now know that there are two distinct 
classes of deformation belts (ridge belts and fracture 
belts) in Lavinia. A particularly important question is 
whether or the not ridge belts and fracture belts both 
owe their origin to a similar mechanism. Evidence sug- 
gesting that they do includes geometric similarities be- 
tween the two classes of belts (widths, positive topo- 
graphic signatures, patterns traced across the plains) 
and, particularly, the observation that a few belts un- 
dergo a transformation from ridge morphology to frac- 
ture morphology along their length. Similar changes in 
morphology along the length of a belt occur in the de- 
formation belts of Atalanta Planitia as well [Solomon et 
al t this issue]. It is clear, however, that whatever the 
underlying cause of belt formation, the surface tectonic 
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expressions of ridge and fracture belts are markedly dif- 
ferent. 

Because of the gross similarities of both types of 
belts, and particularly because both are elevated, one 
reasonable working hypothesis is that both result from 
belt-normal crustal shortening and thickening [Solomon 
ct a/., 1991]. In the case of the ridge belts, this argu- 
ment is straightforward: the ridges are interpreted to be 
folds formed by ridge-normal compression. In the case 
of the fracture belts, however, the argument is really 
based only on the altimetry. All of the belts lie higher 
than the surrounding plains and hence are attributed 
in this interpretation to a thickening, via belt-normal 
shortening, of the buoyant crust. 

A difficulty with the idea that both types of belts 
share a common origin is that the most straightforward 
interpretation of the observed tectonic features in the 
fracture belts runs directly counter to this interprets 
tion. While near-surface materials in the ridge belts 
have been deformed to produce folds, the same mate- 
rials in fracture belts are intensely faulted in what ap- 
pears to be a manifestation of near-surface extension. 
One could interpret these belts, then, as being wholly 
extensional features, with the uplift a consequence of 
lithospheric heating or igneous intrusion. Alternatively, 
in the compressional interpretation, the observed fault- 
ing could be the consequence of stretching of brittle 
surface materials across the crest of an archlike uplift 
whose underlying origin is compressional. The strongest 
argument that both belt types owe their origin to the 
same process may be that single belts exist that exhibit 
both morphologies, depending on their local trend. 

While the great extent and uniformity of the orthog- 
onal tectonic pattern of Lavinia indicate that the near- 
surface materials there have been subjected to a stress 
field of significant spatial and temporal regularity, the 
deformation belts of Lavinia, in contrast, are geometri- 
cally very irregular and can change trend substantially 
along their length, bifurcating and merging to form an 
intricate pattern. If there is a single underlying belt 
formation process, then the irregular orientations of 
the belts with respect to the orthogonal pattern of the 
plains suggest that this process is not related in a sim- 
ple manner to the one responsible for the near-surface 
stresses expressed on the plains. The morphologic mani- 
festation of each belt nonetheless commonly depends on 
the orientation of the belt with respect to the tectonic 
pattern of the plains: ridge belts form where the belt 
is perpendicular to the regional compressive stress, and 
fracture belts form where the belt is perpendicular to 
the regional extensional stress. In both cases there ap- 
parently can be some distributed shear across the belt, 
as indeed is likely to be the case when the belt axis is 
not perfectly orthogonal to any principal stress. 

The source of stress responsible for formation of the 
deformation belts of Lavinia is not known, but man- 
tle convection must be considered the prime candidate. 

It has been suggested that the Venusian mantle may 
lack an Earth-like asthenosphere [e.g., Phillips, 1990]. If 
this is the case, convective stresses can couple strongly 
to the lithosphere and can lead to deformation of the 
lithosphere at wavelengths comparable to the spacing 
of the Lavinia deformation belts [Zuber, 1987, 1990). 


The surface manifestation of convective stresses, how- 
ever, will also depend on the rheologic structure of the 
crust. A model in common use for crustal and upper 
mantle rheology on Venus involves a crust of basaltic 
composition overlying an olivine-rich mantle. This con- 
figuration can lead to a layered rheologic structure, with 
a strong upper crust, a ductile lower crust, a strong up- 
per mantle, and ductile mantle material below [Zubtr, 
1987; Banerdt ond Golombtk , 1988]. A weak lower crust 
may be relatively ineffective in transmitting convective 
stresses to the upper crust. 

We cannot presently answer the question of whether 
the mantle motions and resultant stresses that led to 
formation of the two classes of belts in the Lavinia re- 
gion were fundamentally the same. In some respects 
the interpretation of the data is more straightforward if 
they were not, since the surface morphologic expressions 
of the two belt types are quite different. However, we 
have already noted the arguments for the view that both 
belt types share a common underlying origin, and these 
arguments appear fairly persuasive. Because the pat- 
tern of the belts is so complex, stresses produced in the 
strong upper mantle by convection must also have been 
complex, both spatially and temporally. The stresses 
that produced surface deformation on the plains, how- 
ever, must have been spatially and temporally quite 
regular. The surface stress field must therefore have 
resulted from the superposition of (1) stresses gener- 
ated by mantle convection and transmitted across the 
ductile lower crust with (2) stresses in the strong upper 
crust of unknown origin that were much more regular 
in space and time. The former stresses are expressed at 
the surface primarily by the elevated belts themselves. 
The latter stresses must have dominated the in-plane 
components of stress in near-surface materials, control- 
ling both the orientation of the plains tectonic features 
and the style of deformation within the belts. A mixed 
stress field of this sort might be enabled by a ductile 
lower crust that incompletely transmits stresses from 
the upper mantle to the upper crust. 

Our observations of the Lavinia region raise a number 
of questions: 

1. What process is responsible for the stress field 
that produced the orthogonal tectonic pattern on the 
plains, both here and elsewhere on Venus? 

2. What process is responsible for the stresses that 
produced the elevated deformation belts, and for the 
concentration of the deformation into discrete belts? 

3. What do the dominant wavelengths of deforma- 
tion, both within and among the belts, indicate about 
the mechanical properties of the venusian crust and 
mantle? 

4. If the two classes of belts share a common origin, 
is there a fully self-consistent stress distribution and 
history that can produce the deformation observed? 

5. If the two classes of belts have fundamentally dif- 
ferent origins, why are they so similar in some respects? 
Why are there belts with both characteristics, and espe- 
cially belts that change style along their length? Why 
are the fracture belts different from rift zones elsewhere 
on the planet? 

In the two sections that follow, we make a first at- 
tempt to address questions 2 and 3 above. The others 
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will be left for consideration in future work. 

Dynamical Models for Bell Formation 

We have argued above that the large-scale morphol- 
ogy of deformation belts, and of the ridge belts in par- 
ticular, may be consistent with the hypothesis that 
these features are products of lithospheric shortening 
and crustal thickening. This was also the conclusion 
reached by most investigators prior to the Magellan mis- 
sion, and previous workers have postulated this shorten- 
ing to be the result of convective downwelling centered 
beneath lowlands regions [Zuber, 1990; Binischadlcr et 
a/., 1990; Phillips et a/., 1991]. In this section we use 
finite element models to investigate quantitatively the 
relationship between mantle convection and crustal de- 
formation, with particular emphasis on the conditions 
favoring the formation of deformation belts. 

The geometry and characteristic spacing of the de- 
formation belts of Atalanta Planitia suggest that these 
features result from lithospheric instabilities induced 
by horizontal extensional or compressional stress [^li- 
ter, 1987]. There are two characteristic scales of de- 
formation associated with these belts: one describing 
the spacing of features (e.g., ridges) within the belts, 
and one describing the spacing of the belts themselves. 


Linear stability analysis indicates that two scales of de- 
formation can arise given sufficient compression or ex- 
tension if the lithosphere is Theologically stratified, i.e., 
if a weak lower crustal layer is sandwiched between a rel- 
atively strong upper mantle and a strong upper crust. 
A stress-dependent rheology is required for the instabil- 
ities to grow. Zuber [1990] suggests that the compres- 
sional stresses required to induce the instabilities may 
be supplied by mantle convection. This hypothesis re- 
quires that the stresses associated with the motions of 
the mantle couple into the overlying crust. The high 
correlation of long-wavelength topography and gravity 
and the large apparent depths of compensation suggest 
that Venus lacks an upper-mantle low-viscosity zone 
[Kiefer et aL, 1986; Bills et a/., 1987; Phillips et a/., 
1991; Smrekar and Phillips, 1991]. If so, then mantle 
convection should indeed couple strongly to the litho- 
sphere, and patterns of mantle flow should have recog- 
nizable signatures in surface topography and deforma- 
tion [Phillips, 1986, 1990]. 

On the basis of analytical models of convection- 
induced crustal flow [e.g., Bindschadler and Parmen- 
tier , 1990; Schmeling and Marquart , 1990], the crust 
may be regarded as having two responses to convective 
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Fig. 17. The finite element mesh used for calculations of convectivdy induced deformation of the crust. The mesh 
is 28 by 80 elements with a highest vertical resolution of 4 km in the upper 40 km. 



13,594 c ^ 

Squyres et al.: Plains Tectonism on Venus 

downwelling The first is a “flexural” response, dur- 
ingwhich the surface and crust-mantle boundary can 
bo thought of as deforming “in phase.” Following this 
mode of deformation, the crust responds by thickening, 
and eventually the surface and the crust-mantle bound- 
ary become anticorrelated. 

Simons et al [1991] have proposed that deforma- 
tion belts like the ones observed in Lavinia Planitia 
can form during the transition from the initial flexural 
response to the later phase of crustal thickening, and 
we expand on this suggestion here. We illustrate this 
proposal by means of idealized finite-element models of 
the dynamic interaction between the crust and man- 
tle. The crust and mantle in these models are treated 
as viscous fluids, and a modified version of the two- 
dimensional Cartesian finite-element program ConMan 
j/finy et al, 1990] is employed. A passive marker chain 
tracks the crust-mantle interface and permits variation 
in the crustal buoyancy as well as specification of dif- 
ferent crustal and mantle rheologies. In addition to the 
flow field, the stress field ip the lithosphere and the sur- 
face topography are readily calculated. 

We use an irregular finite element mesh that is 28 
elements high and 80 elements wide (Figure 17). Our 
maximum resolution is in the 40-km-thick top layer 
where each element is 4 km high and 10 km wide. In’ 
aU, the mesh is 800 km in the horizontal dimension and 
400 km m the vertical dimension. We impose free-slip 
boundary conditions on the top and side walls, with 
no flow through these walls. On the bottom boundary 
flow is constrained to be vertical, with no deviatorie 
normal stress permitted. In effect, this last boundary 
condition gives us a virtual 800 km by 800 km box. 

The surface topography is calculated from the vertical 
stresses on the top wall of the box. Top and bottom 
temperatures are fixed at 500*C and 1250 8 C, respec- 
tively. Initially, we impose a linear temperature gradi- 
ent across the thermal boundary layer and set the rest 
of the mantle to be nearly isothermal with a small si- 
nusoidal temperature perturbation imposed throughout 
the box to initiate flow. 

Our exploratory models have a constant viscosity 
mantle overlain by a crust of greater constant viscosity. 

As noted above, the viscosity in the actual Venusian 
crust and mantle is likely to vary with depth; in these 
simple models the higher viscosity crust is intended to 
represent a strong surface layer capable of being de- 
formed by stresses induced by mantle convective mo- 
tions. We adopt a reference mantle viscosity of 10 21 
Pa s, a crustal density of 3000 kg m" 3 , and a mantle 
density of 3300 kg m" 3 . We use initial crustal thick- 
nesses of 12 and 25 km, and raticn of crustal viscosity 
to mantle viscosity of 25, 50, and 100. In all these mod- 
. ' i convection produces horizontal compressive stresses 
in the lithosphere that overlies downwelling mantle and 
extensional stresses in the lithosphere that overlies up- 
welhng mantle. As the convective vigor increases, so 
does the magnitude of the compressive stress over the 
downwelling, with stress in the crust reaching values 
in excess of 100 MPa in about 100 m.y. (Figure 18). 

Note that the rate of increase in compressive stress de- 
creases with increasing crustal viscosity. This is because 
the stronger the crust, the more the development of the 



Fig. 18. Development of horizontal comprewive stress versus time 
in model* of convectiyely induced enutal deformation. Strew i* 
mewured at the top right-hand (downwelling) aide of the model. 
Models shown indude crust/mantle viscosity ratios of 25 (solid 
hne), 50 (dashed line), and 100 (dotted line), and a crustal thick- 
ness h of 12 or 25 km. Note the similarities among the curves 
and the peak in compreasional stress. 


convective instability in the mantle driving the deforma- 
tion is impeded. We also find that the magnitude of the 
peak compressive stress achieved above the downwelling 
increases with higher viscosities and/or with thinner ini- 
tial crustal layers; the stronger the crustal lid, the larger 
the tractions from mantle convection supported by the 
crust. Since force balance on the crust requires that 
shear traction integrated along the base be balanced 
by normal stresses integrated through its thickness, the 
thinner the crust , the larger the horizontal stresses. 

The topography over the downwelling undergoes a 
complex evolution in these models (Figure 19). Initially, 
the effect of stresses induced by downgoing mantle dom- 
inates, producing topographic lows. After reaching a 
minimum, the topography grows due to the isostatic 
contribution from the thickened crust. We suggest that 
if the deformation belts are the result of compressive 
stress and that mantle convection is the source of this 
stress, then the belts will develop at or before the time of 
maximum horizontal compressive stress. This time cor- 
responds to the transition from decreasing to increasing 
topography. 6 

Theory derived by Biot [1961] shows that horizontal 
compression of a buoyant viscous plate overlying a vis- 
cous substratum will have a dominant wavelength that 
depends on the difference in layer densities A p, the ratio 
of viscosities Tjc/rfm, the applied stress P, and the layer 
thickness H c . For values of the quantity P/(ApgH e ) 
less than about 7, the amplification of folds becomes 
insignificant. For one of our models described above, 

P - ?i#MP '??■ = 12 km ' = 300 k«/nA -J 

r - iuu Mr», the dominant wavelength is about 100 
km However, these parameters are near the threshold 
conditions for significant fold amplitudes. 

Our initial results show that coupling of convective 
motions to the lithosphere will cause significant crustal 
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Fijr. 19 . Evolution of surface elevation (top) and the depth to the crust-mantle interface (bottom) with Ume for 
a model of crustal deformation induced by mantle convection. Note the initial decrease in elevation over the zone 
of downwelling (downward arrow), followed by an eventual increase in elevation. Note also the steady mCTease in 
crustal thickness over the zone of downwelling. The initial crustal thickness, H c , is 12 km, and the ratio of crustal 
viscosity to mantle viscosity is 70. 


deformation. This deformation is complex and highly 
time dependent. Our results are preliminary, and mod- 
eling with more realistic rheologies is in progress. How- 
ever, the simple rheologies assumed here are sufficient 
to generate stress fields that, together with simple the- 
ory [Biot, 1961], predict the predominant wavelength of 
deformation seen in Lavinia Planitia. This theory may 
also explain why the belts are only generated for a fi- 
nite period of time before they become overwhelmed by 
larger-scale deformation. 

Buckling Within Individual Ridge Belts 

We have argued that the ridge belts observed on 
Lavinia Planitia are fold belts, formed by buckling of a 
strong upper crust in response to ridge-normal compres- 
sion. If this interpretation is correct, then we may use 
observations of ridge geometry, along with a simple me- 
chanical model of buckling, to place some quantitative 


constraints on the structure of Venus’s crust at the tune 
of ridge formation. We model the ridges here as having 
formed due to a buckling instability in an elastic layer 
with a depth-dependent yield strength subjected to a 
compressive end force. Our general approach follows 
that of Banerdt and Golomlek [1988]. We consider the 
strength of the venusian lithosphere to be determined by 
two deformation mechanisms: brittle failure and ductile 
flow [see Brace and Kohlstedt, 1980]. At a given depth, 
the local pressure and temperature determine the yield 
stress for each mechanism; the lithospheric strength is 
taken to be the weaker of the two. 

In our modeling, the critical differential stress for 
brittle failure is taken to be r = fitr n [ Byerlee , 1978], 
where r and <r n are the shear and normal stresses, re- 
spectively, and /i is the coefficient of friction. The criti- 
cal differential stress for ductile flow is derived from the 
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creep equation i = A<r n exp(— QfRT), where i is the 
strain rate, cr is the differential stress, R is the gas con- 
stant, T is the absolute temperature, and Q and n are 
experimentally determined constants. We take values 
for these constants appropriate for dry diabase [Cans- 
tan , 1982], 

An example of a lithospheric strength profile for hor- 
izontal compression calculated from the above models is 
given by the heavy solid line in Figure 20. The strength 
increases linearly with depth near the surface (in the 
brittle regime) due to the increasing lithostatic stress 
and then decreases exponentially (in the ductile regime) 
due to the thermal gradient. For this figure we use a sur- 
face temperature of 720 K, a thermal gradient (assumed 
constant) of 15 K km" 1 , and a strain rate of 10” 15 s" 1 . 
Now, consider how the strength profile of Figure 20 re- 
sponds to an applied stress of of 70 MPa (the dashed 
vertical line in the figure). Only those depths with a 
lithospheric strength exceeding the applied stress (be- 
tween Z t and Zb) will support the entire stress load 
elastically. Depths above Z t will only support part of 
the load elastically (up to the yield stress) and will un- 
dergo brittle failure. Depths below Zb will only support 
part of the load and will undergo ductile creep. 

At the crust-mantle boundary, a change in compo- 
sition can create another rheologic transition. Mantle 
material on Venus has been assumed in most models 
to be olivine-rich in composition and hence to possess 
much higher yield stresses than the crust, leading to an- 
other brittle-ductile failure envelope in the upper man- 
tle [e.g., Zubcr , 1987; Bantrdt and Golombek, 1988]. In 
line with the calculations above, this deeper strong zone 
could be responsible for longer wavelengths of deforms 
tion on Venus. Here we only consider short- wavelength 
folding and take it to be a function only of the crustal 


failure envelope, neglecting the potential influence that 
the coupling between the crust and mantle discussed in 
the previous section could have on the stress field. 

An elastic plate overlying an inviscid fluid will buckle 
if subjected to a sufficiently large compressive end load. 
The critical force required is given by [Turcoitc and 
Schubert, 1982, p. 124] 


„ _ r Eh^pg 1 


1/2 


( 1 ) 


where E , A, and u are the Young’s modulus, thickness, 
and Poisson’s ratio of the plate, respectively, p is the 
fluid density, and g is the gravitational acceleration. 
Given this force, the elastic plate will deform with a 
characteristic wavelength 


A r Eh 3 1 1/4 

~ 2 * L12(1-^WJ 


( 2 ) 


This simple buckling model was adapted by McAdoo 
and Sandwtll [1985] to the more complicated case 
for which the elastic plate possesses a brittle-ductile 
strength envelope as discussed above. They found that 
equations (1) and (2) hold quite well if the plate thick- 
ness is set equal to the thickness of the “elastic core”, 
defined as the region for which the strength of the litho- 
sphere excedes the applied stress (i.e., the thickness of 
the region between Z% and Z% in Figure 20). This ap- 
proach avoids the unreasonably high critical buckling 
stresses found when finite strength envelope is ignored. 
The critical load for this model is determined by equat- 
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20) due to ductile creep. The bottom of a thick crust 
contributes very little strength to the crust, and there- 
fore has comparatively little influence on the folding me- 
chanics. As the crustal thickness decreases, the folding 
wavelength shows more and more of a dependence, with 
the thinnest crust yielding the shortest wavelengths. Of 
course, the development of a second, longer wavelength 
of deformation does depend strongly on crustal thick- 
ness [Zuber, 1987, 1990]. 

The thermal gradient and strain rate both can have a 
significant effect on the folding wavelength. Increasing 
the thermal gradient decreases the wavelength, while 
increasing the strain rate increases the wavelength. Of 
course, as Figure 21c demonstrates, the folding wave- 
length can remain essentially constant as the combina- 
tion of a relatively large strain rate and thermal gradient 
is varied to a relatively small strain rate and thermal 
gradient. This behavior is a consequence of the form 
of the creep law for the ductile deformation. Increasing 
the thermal gradient and decreasing the strain rate both 
lower the ductile yield stress of the crust, thinning the 
lithosphere and decreasing the resulting folding wave- 
length. 


In Figure 21 d, we have set the folding wavelength 
to 8 km, a value representative of ridge spacing in the 
ridge belts of Lavinia Planitia. This allows us to vary 
all three of the key model parameters simultaneously. 
The contours show the thermal gradient required, in 
kelvins per kilometer, to produce folding with a wave- 
length of 8 km, as a function of strain rate and crustal 
thickness. Ultimately, we would like to be able to use 
this plot to constrain the thermal gradient in Lavinia 
at the time of deformation. In order to do this, how- 
ever, we need to know the strain rate involved in ridge 
belt formation. Given the geometry of the belts and 
typical estimates of the age of the Venusian surface, 
a strain rate corresponding to 10% shortening per 100 
m.y. is probably an extreme lower limit; this value is 
8 x 10* 17 s* 1 . Taking this as the strain rate, Figure 
21 d yields a minimum thermal gradient of ~30 K km* 1 . 
This value is somewhat greater than previous estimates 
obtained by similar arguments [ Zuber , 1987, Bancrdi 
and Golomltk , 1988]. The difference is attributable 
to the ability of Magellan to resolve buckling features 
with smaller wavelengths than were inferred from low- 
resolution data. Although Figure 21 d suggests that the 




Fig. 21. Results of buckling calculations for the Venusian crust, (a) Dominant buckling wavelength as a function of 
strain rate and crustal thickness, assuming a thermal gradient of 25 K Ion* 1 . (6) Dominant buckling wavelength 
as a function of crustal thickness and thermal gradient, assuming a strain rate of 10“ 16 §“*. (c) Dominant 
buckling wavelength as a function of strain rate and thermal gradient, for a crustal thickness of 10 km. (d) 
Thermal gradient as a function of strain rate and crustal thickness for a buckling wavelength (8 km) typical of 
the values observed in the ridge belts of Lavinia Planitia. 
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crust&i thickness is not well constrained, the thickness 
of the elastic core is well constrained, with a value near 
about 200 m over the entire parameter space. 

As noted above, these results assume complete de- 
coupling between the crust and mantle. If some degree 
of coupling does exist, as we assumed in the convection 
model above, the buckling wavelengths presented here 
would increase [Banerdt and Golombek, 1988], leading 
to still higher inferred thermal gradients. Conversely, 
the buckling wavelengths would decrease if lithospheric 
strengths are reduced near the brittle-ductile transition 
by 8emibrittle deformation processes [Kirby, 1980]. Fi- 
nally, we note that there is significant uncertainty in 
the choice of creep law parameters for the Venusian 
crust and that changes in a number of parameters cho- 
sen could have a significant effect on the results. 

Conclusions 

The Magellan spacecraft has provided the first de- 
tailed look at the deformation belts of the Venusian 
lowland plains. In some ways, these observations ap- 
pear consistent with pre-Magellan expectations. The 
belts are elevated zones of intense tectonic deformation, 
and plausible arguments may be made that they are the 
consequence of crustal shortening and thickening. The 
observed scales and spacings of the belts may yield in- 
formation about the strength properties of the lower 
crust and upper mantle, and the spacings of tectonic 
features within the belts may carry information about 
the strength properties of the upper crust [ Zuber , 1987, 
1990]. However, there are some surprises as well, the 
most notable one being the fact that there are two dis- 
tinct types of belt morphology. We ultimately must un- 
derstand the reasons for these differing morphologies, 
and how they may be related to the orthogonal pat- 
tern of smaller deformational features that transects the 
plains. Perhaps the most important part of this process 
will be the examination of future Magellan data. These 
data for Lavinia will include additional SAR images at 
incidence angles that will allow stereographic determi- 
nation of topography, along with high-resolution map- 
ping of the planet’s gravity field. Moreover, images of 
other deformation belts, particularly those in Atalanta 
Planitia, have only recently become available for de- 
tailed study. A systematic comparison of the deforma- 
tion belts throughout Venus will provide a foundation 
for further modeling of belt formation and evolution. 
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